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PREFACE 


The subject matter of this research 
program includes the proper description and 
interpretation of field ionization data, the 
improvement of ionizer designs, and the adapta- 
tions to mass spectrometry. The program is 
sponsored by the National Aeronautics and Space 
Administration, Electronics Research Center, 
Boston, Massachusetts. Dr. Wayne L. Lees is 
the project monitor for NASA-ERC. Mr. Paul R. 
Yeager of the NASA Langley Research Center is 
also providing monitoring services. 

This interim scientific report covers 
the work performed in the first year of the 
program, 25 March 1968 through 25 March 1969. 

The work- has been performed in the Physics 
Department, University of Missouri - Kansas City 
Dr. Paul J. Bryant and Dr. James M. Phillips 
serve as principal investigators and authors of 
this report. They were assisted during this 
report period- by Mr. Paul L. Gutshall, Mr. G. 
David .Moore, and Mr. Charles Gosselin, staff 
members; Mr. Dale Holt, Miss Carol Stephenson 
and Mr. George Mansell, graduate research 
assistants . 

Communications and suggestions regarding 
this research program will be sincerely appre- 
ciated. - 
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SUMMARY 


This report describes the development of an efficient 
field ionization ion source and quadrupole focusing system 
for the production of ions for mass spectrometry. The 
phenomena of field ionization by the quantum mechanical 
" tunneling 1 ' of an electron out of the potential barrier of 
an atom or molecule has been discussed and utilized by a 
number of workers . 1 , 2 , 3 , 6 , 7-/ The application <pf field 
ionization sources to mass spectrometry has been limited 
by the low ion yield as compared to other standard ioniza- 
t ion techn ique s . 

A new carbon filament ionizer, designed and utilized 
on this program, has increased ion* yield by two orders of 
magnitude, making this technique competitive with electron 
impact sources. The new ionizer coupled with an ion lens 
system designed for focusing and retarding the high energy 
field ions has been successfully used as an ion source for 
a commercial quadrupole mass spectrometer. 

The mass spectrometer resulting from the union of the 
commercial unit with the field ion source and lens system 
has several unique features. The gas sampling and spectro- 
meter chambers are separately evacuated so .that toxic and 
corrosive gaseous environments may be analysed without 
damaging or contaminating the mass spectrometer. Gases in 
the sampling chamber may be analyzed by either electron 
bombardment or field ionization and the results compared. 

The field ionization source and lens system has been specifi- 
cally designed to reduce the possibility of a background 
signal due to spurious radiation such as x-rays or electrons. 
The sets of quadrupole lenses and immersion type lenses 
accomplish focusing over- a Tong path length while reducing 
background affects from the ion source - to negligible levels. 

Thus, the high ion yield obtained with the new type 
field ion sources, and the preliminary results with, the mass 
spectrometer show the feasibility of practical applications. 

In addition to the experimental results, ion yields have 
been computed from the transmission probabilities for an 
electron to pass through a molecular potential barrier. 
Transmission -probabilities were calculated by two methods, 
the WKB method, and a direct solution of the linear Schro- 
dinger equation. Both of these methods have been programed 
for the computer so that a direct comparison between experiment 
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and theory is possible by inserting the experimental para- 
meters into the computer programs. The results of a third 
calculation of the transmission probability for the inter- 
mediate field range predicts an increase of ion yield of two 
• orders of magnitude for a gain in field strength of 1 V/A. 

The results .of these calculations will serve as a guide 
to optimization of ion yield in a working mass spectrometer. 
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I. INTRODUCTION 


The build-up of toxic trace contaminants in a cabin 
atmosphere poses a serious problem for long term opera- 
tion. The use of a mass spectrometer employing an elec- 
tron impact ionization source for monitoring the build up 
of these trace contaminants has several disadvantages: 1) 

the cracking patterns complicate the analysis of the 
atmosphere; 2) in many cases parent molecules are eliminated 
by electron impact-; and 3) the telemetering and monitoring 
of such data is complicated by the large number of signals 
and the multiple contributions to individual signals. 

The use of a field ion source for mass spectrometry 
provides solutions to the problems presented above. The 
electron tunneling process in field ionization may occur 
with the least' amount of energy input to a molecule ; thus 
fragmentation and the resulting complexity of the mass 
spectra are minimized. The parent molecule generally yields 
the only significant, or at least the most intense peak 
of the mass spectra. The reduction in "cracking" and the 
intensity of the parent peaks minimizes the problems of 
telemetering and monitoring mass spectrometer data. 

The low yield of typical field ion sources (tips and 
wires), as compared to electron impact sources, has limited 
their usefulness in mass spectrometry. The purpose of the 
research to be described in this report is the development 
of a field ionization source with improved ion yield for 
application to mass spectrometry. 

The purpose of the theoretical calculations is to provide 
support of the experimental effort and to develop' a more • 
fundamental description of the basic physical mechanism of 
field ionization. These goals are -achieved by a study of 
the kinetics of the molecular interactions with the emitter 
geometry and a quantum mechanical theory of the ionization 
rates. The methods employed have included: the WKB 

Solution of Schrbdinger 1 s equation for bound electrons in 
the potential well of a molecule in a high electric field; 
and band, theory considerations for the ionizer material. 

II. EQUIPMENT AND RESULTS 

The experimental portion of the program is presented 
in two sections; (A) ion yield measurements and (B) mass 
spectrometry. The experimental procedures and results for 
the ion yield measurements are presented in. Section A. The 
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design, construction, and testing of special equipment 
necessary for the adaptation of the field ion source to 
mass spectrometry is discussed in Section B. 

A. ION YIELD 


The total ion production from several experimental 
ionizers was measured with special equipment (see section 
1) under known conditions. The data is presented in section 
2 below in comparison with other ionizer designs. 

1 . Equipment 

The equipment required for the ion yield measurements 
includes: (a) a field ion microscope tube, (b) a gas con- 

trol system, and (c) an evacuation system. The overall system 
used for these 'measurements , along with a quadrupole mass 
spectrometer 5 is shown in Fig. 1. 

a. Field Ionization Tube 


A special field ion microscope tube was designed and 
constructed for the measurement of total ion yield. Total 
ion yield values for given- values of gas supply function 
and field strength are necessary for the interpretation of 
field ion data. Thus a spherical tube with a collector 
which completely surrounds the ionizer structure was employed 
(see Fig. 1). 

The tube may be readily disassembled for mounting 
various experimental ionizer structures. Electrical vacuum 
feedthroughs are provided for the high voltage ionizer, 
and the ion collector. An electrical guard ring is also 
provided to prevent current flow from the collector structure 
to the high voltage electrode. A re-entrant glass feed- 
through provides a long path for up to 30 KV insulation. An 
optically transparent stannous oxide coating on the inner 
glass bulb wall constitutes the ion collector. A Keithley 
model 610B electrometer is employed for ion current measure- 
ments. 

b . Gas Control System 

The gas control system consists of the gas inlet unit 
and the pressure measurement system. 

The gas inlet system employs a variable micrometer leak 
valve, connected with stainless steel tubing to a. commercial 
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Fig. 1. Experimental apparatus used to measure total ion yield from 
various field ionizers. 
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vacuum flange. High purity gases from lecture bottles are 
admitted through the variable leak valve. A dynamic equilib- 
rium is established within the system by adjusting the 
micrometer valve until the desired pressure is maintained 
between the inlet flow and the evacuation rate. 


The second part of the gas control system is a dual 
range pressure measurement unit. A Granville-Phillips 
ionization and thermionic gauge control unit is employed 
in conjunction with a Varian model UHV-12 gauge tube. 
Pressure measurements from 10 Torr to lO"^-- 1 - Torr may be 
made with corrections for x-ray photo-current at the lower 
values. Accurate pressure values are needed to determine 
the gas molecule supply function for field ionization. 


c . Evacuation System 


Several different pumping mechanisms are employed for 
various stages of the evacuation process. A mechanical 
fore pump and fluid diffusion pump are used for initial 
evacuation and during fhe period of bake-out of the re- 
maining system. Getter-ion and titanium sublimation pumps 
are available for ultra-high vacuum operation. 

The diffusion pump unit is necessary when an inert gas 
such as helium is being used as the working fluid since the 
getter- ion pump is not capable of handling large quanities 
of inert gases. A special chemical trap is employed to pre- 
vent diffusion pump fluid from back- streaming into the 
system. A new artificial fluid (Santovac 5) compound is being 
used. A large volume titanium sublimation pump is also 
included in the system. 

A differential evacuation capability has been pro- 
vided so that the gas ionization chamber may be exposed 
to some environment of interest while the main body of the 
mas’s spectrometer is maintained in ultra high vacuum. The 
diffusion pump system may evacuate the ionization chamber 
while the getter-ion-sublimation pump system operates on 
the mass spectrometer chamber. Thus, many different gas 
environments including inert and corrosive gases may be 
analyzed ‘while the mass spectrometer is maintained clean 
under high vacuum. 


2. Field Ionization Yield Data 


The total ion current yield characteristics have- been 
studied for experimental field ionizers. Ion current versus 
voltage -data from four different ionizers are presented in 
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Table I. The ions were formed from trace amdunts of the 
toxic gas carbon monoxide. Comparative data for equiva- 
lent conditions is also shown in Table I. The data nor- 
malized to unit gas pressure and unit length of the ex- 
tended ionizers. That is, a comparison of data is pre- 
sented for the number of gas ions formed from an equiva- 
lent number of neutral molecules and from an equal ionizer 
length. Thus the data of Table I may be used to determine 
the actual ion current in amps which can be provided to 
any instrument by multiplying the appropriate ion accept- 
ance aperature dimension and transmission value, and the 
gas pressure by the yield value given in amps/Torr-mm . 

The ion yield values obtained from four of the experi- 
mental ionizers tested under this research program are 
listed under the heading, "filament." The highest yield 
values found in the scientific literature are reported by 
Robertson and Viney2/ using a 10 micron diameter platinum 
wire with CO gas and are listed for various voltages under 
the heading, "wire." The ion yield values obtained under 
this research program with approximately 1 micron diameter 
carbon filaments are from 10^ to 10 3 times higher than the 
best values reported for a wire. In addition, the wire 
data is about 10 2 times greater than values reported for a 
metal tip. Thus the ion yield obtained with the experimen- 
tal filament ionizers developed under this program is suffic- 
ient for practical applications to the operation of instru- 
ments such as mass spectrometers. 

The experimental equipment employed for the ion yield 
measurements is shown in Fig. 1. Trace quantities of carbon 
monoxide were admitted from a lecture bottle, as seen in Fig. 

1. The spherical field ion microscope bulb collects the ion 
current generated by the experimental ionizers. The electro- 
meter shown is connected by a short length of co-axial cable to 
the field ion bulb. The bulb is shown exposed in F ig . 1 but 
during operation it is enclosed in a Faraday cage for electro- 
static shielding. With this experimental apparatus an ionizer 
filament and holder assembly has been developed for use in 
a mass spectrometer. 


Expressions relating the ion current yield to the elec- 
tric field, emitter geometry, emitter and gas temperature, 
emitter material and the species of gas .molecules are pre- 
sented below for both the high and the low field cases. The 
details of the derivation of these equations are presented 
in Appendix I (pp. 26-31). 
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Voltage 


8000 


13000 


16000 


Ion Current 
Wire 17 


8.3 « 10 -8 


5.2 x 10" 7 


Table I. Response data for the carbon filamei 
platinum wires of Ref. 1; the data : 
unit length of the ionizer. 
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i i 

Yield (amj}/ torr - mm) 


Filament 


( 1 ) 

1.4 x 10 -5 

( 2 ) 

5.3 x 10“ 5 

( 3 ) 

4.7 x 1 ( T 6 

( 4 ) 

4.2 x 10 -5 


( 1 ) 

9.2 x 10" 5 

( 2 ) 

1.4 x 10 -4 

( 4 ) ' 

2.0 x 10 -4 

(1 ) 

1.2 x 10 “ 4 

( 4 ) 

3.8 x 10“ 4 


s of this program compared to the 
normalized for gas pressure and • 








For the high field case the total ion current is 
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Figure 2 shows the relative number of ions formed at a • 
distance X from the emitter surface for low, intermediate 
and high fields. 

The total ion current is dependent on the transmission 
probabilities of electrons through the potential barrier of 
the molecule. Transmission probabilities were calculated 
by both the WKB method and a direct solution to the Schrodinger 
equation. The .details of these calculations are presented in 
Appendix II (pp. 32-37). 

The potential function used in the WKB approximation is 
presented in Figure 3. The transmission probabilities cal- 
culated by the WKB method are presented graphically in 
Figure 4 and in Table II as a function of the energy deficit 

of the ion, E.. 

d 

The potential function used for the direct solution of the 
Schrodinger equation is presented- in Figure 5. The transmission 
probabilities as a. function of the ion energy deficit are pre- 
sented in Figure 6 and in tabular form in Table III. Figure 
7 and Table IV indicate the variation of the Transmission 
probability at the critical distance, X , as a function of 
field strength. 

A comparison of the results of the two calculational 
methods shows that the WKB method predicts a transmission 
probability four orders of magnitude higher than the Schro- 
dinger method for the same field strength values. The Schro- 
dinger -method predicts fine structure in the transmission pro- 
bability .which has been observed and explained by Jason.§/ . 

The direct solution of the Schrodinger equation also indicates 
a change in the transmission probability of two orders of 
magnitude when the field strength, in the intermediate field, 
is varied by 1 V/$. 

Proposed improvements to the theory of field ionization, 
utilizing molecular orbital characterization and Band theory, 
are presented in Appendix III (pp. 38-50). 
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X 


d 

Fig. 3. A diagram of the potential energy of an electron in the presence of a 
molecule, a metal surface, and an electric field.. This potential was 
used in the WKE calculation. 
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TABLE 13 


TRANSMISSION PROBABILITY BY THE V/KB METHOD 


in eV 



) 


6.101 

7.250 

8.421 

9.605 

10.80 

12.00 

13.21 
14.43 
15.65 
16.87 
18.10 
19.33 
20.56 
21.79 
23.03 
24.26 
25.50 
26.74 
27.97 

29.22 
30.46 
31.70 
32.94 
34.18 
35.42 
36.67 
37.91 
39.16 
40.40 


4.1 17 x 
2 . 1 48 x 
1 .251 x 
8.163 x 
5 .929 x 
4.713 x 
4.015 x 
3.592 x 
3.322 x 
3 . 1 43 x 
3.019x 
2.931 x 
2 . 866 x 
2 .817 x 
2.780 x 
2.751 x 
2.728 x 
2 .709 x 
2.694 x 
2.681 x 
2 .671 x 
2 .662 x 
2 .655 x 
2 .649 x 
2.643 x 
2 .639 x 
2.635 x 
2.631 x 
2.628 x 


'°-i 

io ; 
10 , 

10 -3 

10 _3 

10 -3 

1 °_3 
10 , 

10 -3 

10 -3 
10 , 

10 -3 
10 ; 

10 , 

10 “” 

10 -3 

10 -3 

10 -3 

10 -3 

10 -3 

i° i 
1 °"^ 

^3 
1 ° * 

1°1 

,0 -3 

10 -3 
10 ~ 

IO ”' 5 
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Fig. 5. The potential energy plot used in the direct solution of the linear 
Schrodinger equation. 




5 10 15 20 25 30 , 35 40 

E d (eV) 

Fig. 6. 'The transmission probability as calculated by the 

■direct solution of the linear Schrodinger equation. 
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TABLE III 


TRANSMISSION PROBABILITY BY THE DIRECT 
SOLUTION OF THE SCHRODINGER EQUATION 


Ed 


in eV 

6.4 

8.0 

10.0 

12.0 

14.0 

16.0 
18.0 
20.0 
22.0 

24.0 

26.0 
28.0 

30.0 

32.0 

34.0 

36.0 

38.0 

40.0 


'll E d ) 

3.69 x 10 
8.77 x 10 
2.03 x 10 
1 . 15 x 10 
4.72 x 10 
6.98 x 10 
1 .25 x 10 
5.44x 10 

8.62 x 10 
6.27 x 10 
7.09 x 10 
5.86 x 10 
7.35 x 10 
5.06 x 10 
7.12 x 10 
5.50 x 10 
5.00 x 10 
6.38 x 10 


-6 

-7 

-7 

-7 

-8 

-8 

-7 

-8 

-8 

-8 

-8 

-8 

-8 

-8 

-8 

-8 

-8 

-8 
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TABLE IV 


TRANSMISSION PROBABILITY AT .THE 


CRITICAL DISTANCE 
Pm V/A l(E dc ) 


2.0 

1 .91 x 10 

2.2 

5.70 x 10 

2.4 

2.08 x 10 

2.6 

6.24 x 10 

2.8 

1 .59 x 10 

3.0 

3.60 x 10 

3.2 

9.14 x 10 

3.4 

1 .76 x 10 

3.6 

3.14 x 10 

3.8 

5.28 x 10 

4.0 

8.42 x 10 

4.2 

1 .29 x 10 

4.4 

1 .89 x 10 

4.6 

2.68 x 10 

4.8 

3.70 x 10 

5.0 

4.97 x 10 


-7 

-7 

-6 

-6 

-5 

-5 

-5 

-4 

-4 

-4 

-4 

-3 

-3 

-3 

-3 

-3 
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B . Mass Spectrometry 


The special equipment: required for the application of 
field ionization to mass spectrometry includes: (1) a field 

ion source and ion focusing system, ( 2 ) a- mass spectrometer, 

(3) an evironmental sampling system, and (4) an ultra-high 
vacuum system. Each of these systems is described in detail 
below. The overall design of the system is represented schemati- 
cally in Fig. 8. The field ion source and ion focusing systems 
are shown in detail in Figs. 8-10. The operation of the total 
system yielding a field ion mass spectra, is shown in Fig. 11. 

1 . Field Ion Source and Ion Focusing System 

The purpose of this phase of the program was the genera- 
tion of ions from a neutral gas and the introduction of those 
ions into a mass spectrometer. A small diameter wire or hemi- 
spherical tip may be used alternately.?/ as an ionizer in this 
design. An array of quadrupole lenses and immersion lenses . 3,4,5 / 
see Fig. 9, focus the 'ions into a quadrupole mass spectrometer. 

The ionizer system is capable of operating with a poten- 
tial difference as high as 40 kilo-volts. Thus, the desired 
field strength values covering low, intermediate, and high 
fields, from 1.5 to 3.5 volts/ Angstrom, may be achieved with a 
wide range of emitter sizes. A high potential gradient is 
maintained between the ionizer and the first acceleration aper- 
ture. The geometry is arranged so that a large fraction of the 
field lines from the ionizer are directed toward this first 
aperture. Thus a high percentage of the ions should pass 
through the large 2.54 cm diameter acceleration aperture. 

Ions generated in the high electric field gain from 10 to 
30 KeV of kinetic energy in passing. from the ionizer to the 
potential of the accelerating aperture. (Their energy may be 
retarded to a few KeV by an intermediate lens.)' X-rays and 
other spurious radiation may be generated. Thus three pro- 
blems must be handled: focusing of high energy ions; retarda- 

tion of their energy before mass spectrometer analysis; and 
reduction of the incoherent background. The system of high 
voltage double quadrupole and immersion type lenses shown in 
Figs. 8 and 9 is designed to focus the high energy ions onto a 
small exit aperture. Since this aperture subtends an angle of 
only 1/3 of a degree from the source, incoherent background 
signals are substantially reduced. That is, the combination 
of a long path length and a small final aperture reduces the 
background of spurious- signals reaching the mass spectrometer. 

The final stage of the lens system consists of a field 
retardation unit designed to minimize the natural defocusing 
effect. The final energy of the ions must be reduced to less 
than 100 eV for proper analysis by the quadrupole mass' 
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Fig. 8. Ion lens assembly, from left to right,: gas inlet, field ionizer and ion 

focusing sub-systems to supply field generated ions for mass spectrometry. 
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Fig. 9 Interior view of vacuum system showing the field ion lens, left 
and quadrupole mass spectrometer, right. 
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Fig. 10. A gas inlet system and field ionizer for 

direct sampling of high pressure atmospheres. 
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spectrometer. The retarding lens system employs a 3.2 mm 
th,ick disc with a 4 0 mesh per cm nickel screen on both faces. 

This grid screen has a transmission value for ions of 68% 
thus permitting a significant fraction of the ions to pass 
through the retarding plane. 

The retarding field acts as a repeller for ions approach- 
ing the retarding plane thus a defocusing of the beam occurs. 

An additional electrode has been placed in front of the retard- 
ing lens to provide a focusing field for the decelerating ions. 
Following the retarding plane a unipotential lens is used to 
focus the ions through a set of x and y deflection electrodes 
onto the entrance aperture of the quadrupole mass spectro- 
meter. 

2 . The Mass Spectrometer 

Ions - produced by means of "high" or "intermediate" fields 
(from- 2 x 10 8 to 5 x 10 8 v/cra) are found to possess an energy 
Spread of approximately 30 eV . Ions generated by means of a 
"Low" field (1 - 2 x 10 8 v/cm) may be essentially mono -energetic . 
An ionization source for mass spectrometry of gas molecules 
generally requires the use of intermediate or high fields to 
avoid the complications which occur from gas molecule inter- 
actions with the ionizer surface at low field. Therefore, 
the spectrum of ion energies must be tolerated if the complica- 
tions inherent in low field operation are to be avoided. This 
energy spread will not preclude good resolution by a quadrupole 
mass filter type spectrometer which is relatively insensitive 
to energy. 

The spectrometer selected for this study filters the ions 
according to their mass rather than their momentum. This fea- 
ture is particularly applicable for the analysis of field 
ionized molecules which have a fixed mass but varying momen- 
tum values. The original design of such an instrument was 
called a mass filter. That is, an instrument which filters 
selected ions in the region between four equi-spaced rods or 
poles. Thus, the term quadrupole was adopted to describe such 
instruments. 

The commercial unit, purchased by the University of . 
Missouri for this program, is a Spectra Scan 750 constructed 
by Finnigan Instruments and distributed by the Granville- 
Phillips Co. The electronic control unit and the quadrupole 
sensor are shown attached to the vacuum system in Fig.- 1. This 
instrument has a mass- range from 0 to 75 0 atomic mass units 
and will accomplish unit resolution up to 500 amu. Thus the 
toxic spacecraft contaminants of interest in this program lie 
within the working range of this spectrometer. 
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3 . Environmental Sampling System 


The gas inlet and sampling system shown in Fig. 10 was 
designed for direct analysis of gaseous environments over a 
wide range of pressures. Another feature of this design is 
the ability to utilize a high pressure in the ion source 
chamber while the mass spectrometer is under ultra-high 
vacuum (see Section 4 below for a detailed description). ' The 
ultimate advantage of this procedure is an overall increase 
in sensitivity since the ion source operates at high effic- 
iency for the ions because the detector unit is free of inter- 
fering gas. The response of the spectrometer is also stable 
over a period of time and thus available for continuous or re- 
peated comparisons of data. 

4 . Ultra-High Vacuum System 

Two vacuum systems are employed with this mass, spectro- 
meter system.' The first is a high, throughput diffusion pump 
system attached directly to the ion source section. This 
pumping system permits the sampling of any gas, even inert 
species, since the diffusion pump is capable of handling all 
gases. A liquid nitrogen cooled chevron type baffle is used 
to block the backstreaming of hydrocarbons from the pump. 
Another precaution has been taken to protect the cleanliness 
of the ion source " system; a large conductance valve has been 
located between the pumping system and the ion source chamber. 
Thus the chamber may be isolated from the dif fusion_pump dur- 
ing hold periods when the liquid nitrogen trap is not active. 
During these periods of time the chamber may be connected, 
via a second valve,- to an ultra-high vacuum getter-ion and 
sublimation pumping system, see Fig. 7. The latter pumps pro- 
vide continuous ultra-high evacuation of the mass spectrometer 
chamber. When the latter of the two valves, mentioned above, 
is closed the ion spurce chamber is essentially isolated 
from the ultra-high vacuum mass spectrometer chamber and 
opened via the first valve to the diffusion pump. Thus re- 
latively high gas pressures and corrosive gases may be sampled 
without contamination of the mass spectrometer. This ability 
to differentially pump the ion source and detector stages per- 
mits continuous and unrestricted sampling of gaseous environ- 
ments. 

5 . Hass Spectra - Results and Discussion 

The results of. this program to date have established the 
feasibility of utilizing a carbon filament as a field ion 
emitter source for mass spectrometry. The improvement in ion 
yield with this new design ion source should' make this tech- 
nique competitive with conventional electron bombardment mass 
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spectrometers . 


The ion lens and energy retarding system are currently 
being perfected. Preliminary mass spectra, see Fig. 11, 
have been obtained with the complete field ionization mass 
spectrometer system. The mass spectrum shown in Fig. 11, 
although preliminary, does show that trace signals (10“^-° amp) 
can be distinguished from the background. A recording sys- 
tem biased above the background level (indicated by the 
dashed line in Fig. 11) would respond to the two trace gas 
signals for mass numbers 18 and 28 as seen in Fig. 11. 


Total ion yield measurements in the bulb will be con- 
tinued in an effort to optimize operating conditions for 
maximum yield and correlate theoretical - and experimental 
results . 
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Fig, 11. Field Ion mass spectra showing trace gas signals 
for water vapor (l8) and nitrogen (28) above 
the backbround level (dashed line). 
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AEPENDIX I. Field Ionization Using a Cylindrical Emitter 

The physical process involved in the generation of ion 
current by field ionization is a complex interaction of 
several variables. The primary variables treated here are: 
electric field, emitter geometry, emitter and gas tempera- . 
rures, emitter material, and the specific gas molecules 
being ionized. 

At high field strengths all of the molecules arriving 
at the emitter are ionized. The ion current, for the high field 
case, is therefore just the ionic charge times the number of 
molecules arriving at the emitter. The number of molecules 
attracted to the emitter is given by the supply function Z 
which exceeds the normal number bombarding the geometric 
area of the emitter because neutral molecules are polarized 
and attracted to the emitter by the high field. At lower 
fields there are many other considerations. The high field 
case is important because higher fields produce greater ion 
currents and the supply function may be checked directly by 
experiment without recourse to transmission probabilities. 

A. The Electric Field of a Cylindrical Emitter 

Previous report have presented analyses, of the 
electric field on the basis of tip emitters. Now we wish 
to consider an emitter with cylindrical geometry. The field 
will cause a polarization of the gas molecules and draw 
them toward the emitter and the region of high ionization 
probability. The emitter is considered to be a cylinder of 
length L and very small radius R 0 with supports at each end. 

The surface of the emitter is considered to' be a "smooth" 
cylindrical surface. The field is assumed to be that of an 
infinite conducting cylinder of radius R c \ The finite length 
of the cylinder and. the end supports are ignored. The dis- 
tances that are important in the capture of molecules are quite 
close to the emitter so end effects contribute little to the 
field over most of the emitting surface. In addition, the 
supports will not have sufficient field strength to contri- 
bute to the supply of molecules. Effects due to image 
forces are not considered in the field contributing to the 
supply function because we are dealing with weakly inter- 
acting forces between the emitter and the polarized molecules. 
However, we cannot neglect these forces when dealing with 
the transmission probabilities of the electrons from the 
molecules to the emitter surface because there the image force 
is between the full charge of the electron and the surface. 
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The field of an infinite conducting cylinder is given 


F = 


0 

£ 



( 1 ) 


where o is the charge in coulombs per unit area. If b is 
the radius of the accelerator (a concentric cylinder to the 
emitter) then at r = b the potential is equal to zero. At 
r = the potential is V yielding 


V, 


R In 
o 


b 

R 


= F r 


( 2 ) 


Previous work with metal 'tips was concerned with an 
image field (for field emission) and the compression factor 
due to the shank effects upon the 'tip end. These items are 
not considered important for cylindrical emitters, for the 
reasons described above. 

E. The Supply Function of Gas Molecules to a Cylindrical 
Emitter 


The potential energy for a particle in a field F is 


-V ( F ) = P ,F + igaF 2 (3) 

d 

where a is the polarizability of the gas molecules and P^ is 
the permanent dipole moment. 


Since the field is a function of r only for this central 
force problem, the angular momentum is conserved. Such a 
problem can be treated as one-dimensional with the centri- 
petal force taken into account. This is done by forming an 
effective potential. The expression for total energy is 

E = T + V = kT = % mr 2 + mr 2 ^ 2 + V (r) (4) 


The angular momentum is 

2 * 

L = mr <j> , 

The effective potential becomes 


V gff (r) = V(r) + 


2mr 


(5) 
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The angular momentum can be expressed as 
L = fflvr 


where r" is the distance of closest approach if the field 
were zero. The term v is given by 


v 


v 

rms 



( 6 ) 


where the factor 2 is due to the cylindrical geometry which 
has two dimensions of escape, while the spherical geometry 
has three. The angular momentum is now 

(2mkT)^r' (7) 


Substituting for L in Eq. 5 we have 

V f f = V(r) + kT. 


( 2 
r 


( 8 ) 


dr 


The centripetal velocity is - 

-r ='V1( E - V eff) % • 


thus from Eg. 4 we have 


For impingement, r must be real at R , and for grazing inci- 
dence at r = R , 
o 5 


E 


V 


ef f 


0 . 


For cylindrical geometry a Z-component of velocity cannot 
aid the escape of a gas molecule from the electric field. Thus 
the Z-component of energy E 3 = 3/2 kT will be taken to be 
E 2 = kT since E„ is parallel to the long cylinder axis. We 
may write 



= 1 


V(r') 

kT 


G 


C 


(9) 


noting that -V(r) is positive. The effective cross-sectional 
area of capture is 2wrL where L is the length of the cylindri- 
cal emitter. The supply function for the cylindrical emitter 
is 


Z = 2irR a hL . (10) 

o c 

If the ionic charge is q and the number of molecules 
hitting a unit cross-section of area per unit time is p(2irmkT) 2 , 
then the ion current for a sufficiently high field is 

i = qZ = q2mR o a c %Lp(2TrmkT) ^ (11) 
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where r is the gas pressure. This Expression (11) is sub- 
ject to the condition that the field is high enough to 
ionize all molecules captured before they reach the cut- 
off distance Xc • Ionization cannot take place for any 
molecule that is closer than X to the emitter surface be- 
cause of the filled states below the Fermi level of the 
emitter material. 


The differences between the low and high field cases 
are illustrated in Fig. 2 which shows the number of mole- 
cules N(X) ionized 'at a distance X from the surface. The 

critical distance X is indicated, 

c 

In the low field case the ionization rate is small 
compared to the rate of arrival of the gas molecules. The 
ion current in this case is the equilibrium concentration 
of gas molecules near the tip divided by the mean life time 
of the ionization of these molecules. 

The equilibrium concentration near the tip differs from 
the concentration .of the gas because of the possible tempera- 
ture difference and the attraction of the polarized mole- 
cules in the high field near the emitter. The tip concen- 
tration can be written as 

C ± - C g f(T g /T t )e“ V( ' F)/kT g (12) 

The function f(Tg/T-j-) is dependent upon the kinetic, 
mechanism. The Boltzmann factor of the field potential is 
analogous to the atmospheric barometer relationship. The - 
function f(Tg/Tp) is is unity at a distance slightly in excess 
of the average hopping height, the temperatures are nearly 
equal at this height. 

Consider a surface at temperature T-t and a gas at tem- 
perature Tg. An equilibrium condition is reached, we now 
seek the zero field equilibrium concentration at the surface. 
Let there be an artificial barrier through which an equal 
number of molecules pass traveling each way. The flux $ 
across the barrier is expressed as 


0 . - $ • 
m out 

and assuming an ideal gas 

P P. 


0 


/ 2 irmkT / 27rmkT , 

g t 
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with 


P = CkT. 


We have 


C 

g 




exp [-V (F )/kT ] 
& 


by adding the effect of the field. 

If full thermal accommodation takes place this relationship 
should hold in the region where most of the ionization takes 
place for our low field restriction. If the concentration 
near the emitter is truly a-n equilibrium concentration, the 
ion current generated in a volume element about the cylindri- 
cal surface is 

di - qC t T -1 dV 

with dV = 2frrLdr and with t the lifetime of an electron in 
its bound state. The calculation of i is the major part of 
the problem and much attention is given to it later. 

We define t as 


T = (v D) 
e 


-1 


(13) 


where is -the arrival rate of the electron to the barrier 
and D i s the tunneling probability. 

For the low field case the total ion current is 


, [T 

= 2TrLC g qx J r f ^ 


R 


exp[-V(r)/kT ]dr 


Letting all functions be at the value they have at R , -as a 
low field assumption, the ion current becomes 


i « 2irLC t _1 qR 


(T 


si* 


R +X 


exp[-V(F )/kT Adr 
o SJ' 

R 

o 
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1 


exp[-V (F )/kT ]X 
o • t c 


(lit) 


2ttLC x _1 qR 
g J o 


It is logical to ask at this point the value of taking data 
at such a low field; some reasons are listed below: 

1) The surface uniformity, i.e. the percent of the 
surface that is actually emitting, when compared 
to the high field case, may be determined; 

2) A measurement when a retarding, potential is used 
of the theoretical X is possible; 

3) The presence of liquid layer may be estimated 
for various emitter temperatures. 

The field strength is high enough in the intermediate 
field case to decrease the steady state concentration near 
the emitter, yet not sufficient to ionize all of the captured 
molecules before they reach the critical distance X . Thus 
kinetic considerations must be used in computing the ion 
current for intermediate field values. Multiple rebounds of 
molecules from the emitter surface with the thermal accommoda- 
tion of each bounce must also be considered . 6 > 7/ Thus, de- 
tailed theories of thermal accommodate on g ? 9,10/ mus t be used 
to predict quantitative values of ion current for the inter- 
mediate field case. 
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AFPENDIX II. Calculation of Transmission Probabilities 

The problem of calculating ion yield currents from a 
field emitter can be divided into two parts: the calcula- 

tion of the supply of ga-s molecules, and the ionization 
rate. The supply of the gas molecules to the emitting' sur- 
face is a problem of gas kinetics and geometry. The 
ionization rate involves the basic quantum mechanical mech- 
anism of the field ion process. 

This report is primarily concerned with the latter pro- 
blem of ionization rate. The ion current from a cylindrical 
emitter for low field approximations is 


i « 2nLC t 1 oR (T /I J^exp[-V(F )/kT JX 
g. -o gt r o g e 


where L is the length of the cylinder, Cg is the concentra- 
tion of gas, t is the life time of the ionization process, 
a is the electronic charge, R Q is the radius of the cylinder, 
Tgis the gas temperature, Tt is the temperature of the 
emitter, V(F 0 ) is the potential about the emitter, the elec- 
tric field is F 0 , and X c is the critical ionization distance. 
We define t as 



where is the arrival rate of the electron to the barrier 
and T is the tunneling probability. 

For the case of a high electric field— the ion current 
is given by 

i =-a2irR o ^Lp( 2irmkT) 22 
o c r 

where 


a = 1 V /kT . 
c o g 

The remaining problem is the solution for the quantum 
mechanical transmission probability. This problem was posed 
in the first report. We restate the problem here with the 
newly obtained computer calculations. The ionization pro- 
bability D is at the heart of this field ionization problem. 
The original linear WKE approach has been applied to the 
problem with many improvement s fi ? 7 , 11 / as described below. 

The condition for the application of the WKB method 
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[hir 3 (dX/dx) << 1 


where dX/dx is the spacial derivative of the de Broglie 
wave length. The above condition is often violated in 
field ionization problems. However, depending on the 
barrier shape, the.WKB method can be valid for a much -j^/ 
wider class of problems than the above condition allows — 


Several three dimensional improvements to the cal- 
culation of field ionization transmission probabilities 
have been given using collision theory!// , variational 
methods!!/ , and perturbation theory!!/. The results pre- 
dicted the half-width of the peak representing- the ion 
energy spectrum. 


Refinements ’ of the image potential have been given — 
as it applies to the reflection and transmission of elec- 
trons through surface barriers. Many of these refinements 
apply to a conceptually similar problem in the periodic de- 
viations from the Schottky line. 


The basic understanding of the interaction between gas 
atoms and metal atoms of the tip in the presence of a high 
electric field took a different turn when structure was 
observed in the energy spectrum of the io n s !° ■ ? — . This 
structure had not been predicted by the quasi-classical WKB 
approximation. Calculations of this structure have been 
undertaken with the following interpretations: the effect 

arises from the fact that the surface potential is resonant 
to the energy of electrons entering the metal;?/ ; and the 
energy spectrum is due to the inference between waves re- 
flected from the rapid variations in the atomic potential 
of the ionizing atom and the waves reflected from ‘the rapid 
change in the potential at the metal surface!// 

The above interpretations of the structure in the energy 
spectrum are, in many ways, quite similar. Conceptual differ 
ences do exist however. These differences in interpretation 
appear to be related in a quantum mechanical sense to the 
Gamow and Fermi theories of alpha decay tunneling. In the 
former theory the tunneling particle was considered as - a 
traveling wave transmitted and reflected from a potential 
barrier while the latter theory considered the particle to be 
initially in a bound state and then tunneling to a free state 

The experimental portion of this program is concerned 
with the improvement of ion current yield. Therefore a 
correlation will be sought with the theory which gives the 
best yield and peak half-width values rather than the best 
description of the structure refinement. 
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1 . Transmission Probability Calculations by the WKB Method 


Comparative calculations have been done to place the var- / 
ious theories into the same format. The linear WKB method^.? 
has been programed for a numerical calculation with a variety 
of conditions. The potential function of Fig. 3 is 


V (x) 




+ Fx- - 


e 

4X 


+ 


e 

£+x 


where x is the distance from the surface to the electron, a 
is the radius, and d is the distance to the ion center. The 
transmission probability is given by 


D(£) = exp 



{V (x) -Fd + Vj 



The results of this calculation are given in Fig. 4 and Table 
II. The calculations were performed on the University of 
Missouri-Kansas City IBM 360 Computer using the Newton- 
Raphson method to find the integration limits and Simpson's 
rule for the integration. In order that these values may be 
compared with the other methods we have chosen the energy 
deficit parameter for comparison 

E d = 0 + P m + Fd - Vi + e/4d 

where 0 =.4.5 eV, p m = 6.4 eV , F = 2.5 eV/A, Vj = 15.6 eV , 
and e is the electronic charge. These conditions are in- 
tended to simulate an experiment with H 2 as the gas and 
tungsten -as the emitter metal. The peak half width (1 eV) 
is in agreement with Tsong and Muller' siiL/ calculated value 
which is twice their experimental value. Boudreaux ^- 5/ • 

using a different potential value and calculation method on 
the one dimensional problem obtained a value four times the 
experimental value.. These authors do not give absolute 
magnitudes for the transmission probability. Jason 3/ does 
give a comparable result, however there seems to be some 
question regarding the units used. In general the WKB cal- 
culation 'is consistent with the usual difficulties, namely 
the half width is twice too wide, no fine structure is 
observed, jJujkz li/ and the magnitudes are far too great. 

B. The Direct Solution of the Linear Schrod'lnger Equation 

A computer program for the calculation of transmission 
probabilities has been written and its limits verified. The 
range of values necessary for directly solving the linear 
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Schrodinger equation has been used for a .potential barrier 
problem. Expressions for the transmission probabilities for 
the barrier in Fig. 5 have been obtained by solving 

d H f 

2-1 - £¥ = o 
d*T 

with 

5 = Tx - (E - $ - yJ/F] 


We have derived the following relationship for the 
transmission probability* which agrees exactly with Alferieff 
and DukeiZ/ 


2mF 




T(E d ) 



2 , 2 -2 . n 2 J 

it P + Q 


where 


P - [Bit5 a )Ai(5 0 ) - Bi(5 0 )Ai(5 d )] 


— [Bi(C a )Ai(e o ) - Bi(5 0 )M(C d )], 


k L k R 


and 

Jr 

Q = ^[Ei(£ d )Ai(5 o ) - Ai(C d )Bi(5 o )] 

R 

+ ^CBi(5 0 )Ai(£ d ) - AiU o )Bi(S d )] . 

L 

where Ai(x) and Ei(x) represent the Airy functions* and the 
prime indicates the derivative with respect to the argument, 
with 


kp = (2mF/fi 2 ) 1/3 s 
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d 

9 

/2m 

V v s 

~ 7 & 


35 



V_5 - 4> + p + F , - 2V T , 

R p m d I 3 


? = - k,-, (E , - 4> - p )/F , 

F ' d T 'm J 


F d 


m ' 


and 


«d * k p [d - <E a - * - % )/F3 


The same parameters were chosen as in the WKB calculation for 
the purpose of comparison. These results are given in Fig. 6 
and Table III. 

The results obtained by this calculation are in complete 
agreement with Alferieff and Duke.ll/ The fine structure 
peaks are positioned near those predicted by the triangular 
well calculation..]/]/ This calculation predicts values some 
four orders of magnitude less than the WKB method. 

The calculations described above will be applied to the 
experimental program. A value of life time (x) will be cal- 
culated from the various experimental parameters. In the 
low field case the height of the main peak at the critical 
distance T(E. ) will give the lifetime. This lifetime with 
the known supply function permits a calculation of the ion 
yield. These example calculations illustrate the workability 
of the method. The computer programs are available and 
working. Only the .proper data cards need be inserted for 
the calculations relationship to the experimental part of the 
program. 

In order to predict the effects of increasing field 
strength on the magnitude of the transmission probability, 
a calculation was performed using the direct solution of the 
Schrodinger equation with a varying field strength. These 
results are given in Fig. T and Table IV. In the range of 
practical experiments there is an increase of two orders of 
magnitude in the transmission probability for each additional 
increase in field strength of .one volt per angstrom. An 
increase of the supply function would increase the ion yield 
even further. 

The calculations reported above have all been c'arried 
out in one dimension. Three dimensional corrections can be 
used to approximate the real situation when required by the 
conditions. This correction can be done by calculating the 
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ionization probability for each 0, neglecting the angular 
Laplacian, then evaluating the integral of the weighted 
probability 


1 

T 



O 

for the total lifetime. The above correction is only 
necessary for the intermediate field case and will not be 
used for low field calculations. 
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AFFENDIX III. Improvements to Field Ionization Theory 

A new theoretical description of the field ionization 
process is now being developed as outlined below. The results 
and predictions obtained should give a better understanding 
of the practical applications of field ionization. 

A. Molecular Orbital Characterization of an Emitter Surface 


A great number of improvements have been offered in the 
theory of the field ionization process and further refine- 
ments would likely prove to be unobservable in the experi- 
mental data. It is striking to notice that all of the com- 
putational models used to date have characterized the emitter 
material by a work function and a Fermi level only. Recently 
published experiments^/ are interpreted to suggest that this 
characterization may not be sufficient. Greater brightness 
in the field ion microscope images of certain crystal planes 
of the emitter suggest that the atomic orbitals that are 
intersected by the surface plane accept electrons from the 
molecules being ionized. The planes which have more orbitals 
protruding ionized more efficiently. In the perspective of 
this model, "field ion images are interpreted as projections 
of regions where the fully occupied orbitals of inert image 
gas atoms can easily overlap with the partially occupied 
single or hybridized orbitals of the surface atoms of the 
emitter . 11 20/ These surface atoms are viewed as preserving, 
to some extent, their individual character. The wave func- 
tions of the surface atoms approximate the wave function the 
same as the isolated atom would have. This view is taken by 
many in the study of chemisorption and catalysis. More support 
for this model is inferred in current theories on hydrogen 
promotion of field ionization and rearrangement of surface 
charge . 21/ 

The unique feature of the experimental phase of the pro- 
gram is the cylindrical emitter of large area and smooth sur- 
face. It would seem natural to use this feature of the emitter 
to study the surface of graphite and improve the theory from 
that direction rather than to concentrate on the field and 
gas molecule side of the problem. 

1. Localized Wave Functions - Delocalized Wave Functions 


Caution must be exercised in considering the molecular 
orbital model of the surface. The electronic configura- 
tion of a metal-adsorbate system in a high electrid field 
has been treated by self-consistent calculations*^/ where 
the possibility of localized surface states are ignored. 
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It is felt that in the field ionization problem the 
separation distance between the gas molecule and the emitter 
surface is sufficiently great so that the proper description 
lies somewhere in between the totally localized and the 
totally delocalized picture of the electronic wave functions. 

2 . The Use of Time Dependent Perturbation Theory in 
Computing Transition Probabilities 

The approach to computing the transition probability 
of an electron from the gas molecule to the emitter surface 
should be attempted by the use of time dependent perturba- 
tion t he or y . 23,15/ 


When the gas molecule is near the emitter surface, the 
electron is subject to perturbations by the surface forces 
and the external field which can induce a transition of the 
electron into the emitter. Letting V be the potential energy 
of all such forces, the probability per unit time of such a 
transition can be calculated using the well-known result 


P 
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g(E e ) 
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U o" 


2 


where $ is the protuding molecular orbiral of the emitter 
surfacef u is the ground sxate wave function of the gas 
molecule , and g(E^) is the density of states. 


The transition time is given by 


The probability of being ionized in the distance L to Z + d 
d s 


PU)dft 


dt _ d£ 

tOl) ~ vt ( i) 


wn ere v is the radial component of the velocity of the mole- 
cule about the emitter*. The probabdfLity of being ionized on 
one pass d.s 


P 

' total 



PU)d£ 


lor a smooth graphite surface the carbon atoms 'are found 
in che usual hexagonal arrangement with hybridized (sp) 2 orbi- 
-ais forming the situna bonds between them. The unhybridized 
parxially filled orbitals are perpendicular to the sur- 
face formd.ng the p.x cloud that rings each hexagon. In this 
approach it is assumed that the transition of the electron of 
she gas molecule will be into this pd cloud. 
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This orbital view calls for a reinterpretation of 
many of the parameters of the field ionization process, 
notably the critical distance x^. A possible inclusion 
of the band structure of the emitter material, graphite , 
may be necessary. 


B. Band Theory and Field Ionization 

The essential physical process of field ionization has 
had several different descriptions . Improve- 

ments to field ion theory have been in two major categories, 
refinements in the potential of the electric field and in 
molecule combination and methods of calculating transmission 
probabilities. Previous theories of field ionization have 
assumed a f ree-electron model for the metallic emitter. The , 
purpose of the present effort is to consider the metal to 
be a system of periodic potentials (see Fig. 12). In such 
a system, the electrons are described by Bloch wave functions 
and the values of real energy are' characterized by a series 
of bands (see Fig. 13). 


The need for a more rigorous theory is indicated by 
recent experimental facts. Knor and Mu 1 1 er.2 0/ observed 
regional brightness in the field ion microscopy patterns. 
Different crystal faces have shown different ionization 
•efficiencies. These results indicate that the periodic 
nature of the metallic crystal is an important consideration. 
An energy analysis of the ion beam by Jason3/ has .shown the 
onset peak to consist of double peaks whose relative heights 
are voltage dependent. Peak height and' resolution are ob- 
served to vary with crystal orientation. 

or an/ 

Surface states (Tamm states ) -? ■ > ■■ result from the 
termination of the periodic system at the metal interface. 
These states are dependent upon crystal direction. They' 
provide energy states for the reception of electrons from the 
gas phase molecules. The transition of electrons from the 
bound stave of the molecule to the Bloch bands of the crystal 
is probably closer to describing the essential physical pro- 
cess than any of the previous theories. 

C. One-Dimensional Brillouin • Zones 


The understanding of a complex quantum mechanical pro- 
cess can be best achieved through a one-dimensional model. 
Once the basic formulation is obtained, ' the extension to 
directional features can be straight forward. An outline of 
the one-dimensional formulation is given in this section. 
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Fig. 13. Energy as a function of wave number for electrons in a one-dimensional 
periodic potential. 


The motion of an electron in a one-dimensional lattice 
is described by the Schrodinger equation with a periodic 
potential (see Fig. 12). We may write it as 

[ - d 2 /dx 2 + V(x ) ] iKx) = Ft#; ( x ) (5) 


where 


V (x + a) = V( x ) . 

The Bloch wave solutions for the allowed values of real 
E are of the form 

^n 3 k (x) = u n 5 k (x)elk?j n - 0, 1, 2, . . . (6) 

where u n ^ (x ) has the same periodicity as V(x) and k is taken 
as real. 1 * The corresponding eigen values are the energy bands., 

E — E^ n - 0, 1, 2, . . . (7) 

The energy bands are separated by forbidden ranges of E 
(see Fig. 13) in which wave functions are .of the form of 
Ec. (6), but with k as a complex variable 

k = - k r + lk i ’ 

When E lies in a band gap there are still two solutions for 
ijj in the periodic potential. An ordinary second order 
differential equation must always have two solutions. These 
solutions are evanescent waves, one increasing and the other 
decreasing on the average exponentially with x. In Parti- 
cular, we have 

t|^(x+a)-= - exp(±k^a)'P + (x) 

for the band gap at k = ir/a 

E may be defined for complex k. Tee energy E is, in 
general complex but it is real along the real axis and at 
values of- k of the form 

k = ( it / a ) ± ik . 

(See Figs. 14 and 15). On the path of real' E one may pass 
continuously from one energy band to another. 
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Fig. 14, Lines of real E in complex }< space. 
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Fig. 15. Regions in complex k space where E is real 
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In consideration of the hulk properties of metals, 
both exponential wave functions are ignored because the 
increasing wave function is not normalized and the de- 
creasing wave function does not extend far enough from 
an individual lattice site to contribute to bulk effects. 

The exponential tails of the metal wave functions are impor- 
tant, (see Fig. 16). At an interface or free surface the 
expontential tails are the extensions of the metal wave 
functions into the vacuum. Since, the exponential tails 
are extensions of the metal wave functions they depend 
upon the k-vector direction, i.e. they vary in length de- 
pending on the crystal direction of the interface. This 
directional dependence is the element needed to explain, 
hopefully, the recent experimental results. Z >12 / 

Fundamentally, the problem is one of matching the 
wave function of the crystal to the wave function outside 
the surface. In the case of field ionization the potential 
outside the surface is further complicated by the presence 
of gas molecules and a high electric field (see Fig. 17). 

The wave ' function for the free side of the surface, is 
for an electron of energy E, determined by the field strength 
and the distance of the gas molecule from the metal surface. 
In the electric field the wave function is in terms of 
Airy functions with E the independent variable. The free 
electron wave function is them matched to the appropriate 
energy Bloch wave function in the metal. 


E and ip can be derived, for complex k, from what is 
known regarding the ordinary band states at real k. In 
fact, the whole modern apparatus of orthogonalized plane 
waves and pseudo-potentials for calculating band struc- 
tures can be applied to determining the evanescent waves . 

With the evanescent waves determined for the periodic cry- 
stal, the matching problem of the \p p (the crystal wave 
function to the i|/j> (the free particle wave function) can be 
carried out. There is only one constant to adjust; there- 
fore, the match is made for a particular E. Plotting 
for and ip]_, over a range of E, an intersection can occur 
thereby determining the allowed surface state (see Fig. 18). 

The wave function inside -the crystal is determined by 
the energy band structure with evanescent solutions of the 
Schrodinger equation. The wave function outside is a sum of 
incident and reflected electron waves. The correct combina- 
tion of these is determined by the matching calculation. 
Electrons incident on a crystal surface are partly trans- 
mitted to states of the same energy inside the crystal and 
partly reflected to states of the same energy outside. The 
knowledge of the correct combination would give the ionization 
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Fig. 16. Potential near a surface. 
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0) In 


ff iciency . Thus, this new theoretical description of the 
ield ionization process should yield data to predict the 
practical applications of field ionization. 
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